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ABSTRACT 

Recent numerical simulations of sunspots suggest that overturning convection is responsible for the existence 
of penumbral filaments and the Evershed flow, but there is little observational evidence of this process. Here 
we carry out a spectroscopic search for small-scale convective motions in the penumbra of a sunspot located 5° 
away from the disk center. The position of the spot is very favorable for the detection of overturning downflows 
at the edges of penumbral filaments. Our analysis is based on measurements of the Fe I 709.0 nm line taken 
with the Littrow spectrograph of the Swedish 1 m Solar Telescope under excellent seeing conditions. We 
compute line bisectors at different intensity levels and derive Doppler velocities from them. The velocities are 
calibrated using a nearby telluric line, with systematic errors smaller than 150 m s -1 . Deep in the photosphere, 
as sampled by the bisectors at the 80%-88% intensity levels, we always observe blueshifts or zero velocities. 
The maximum blueshifts reach 1.2 km s" 1 and tend to be cospatial with bright penumbral filaments. In the 
line core we detect blueshifts for the most part, with small velocities not exceeding 300 m s" 1 . Redshifts also 
occur, but at the level of 100-150 m s" 1 , and only occasionally. The fact that they are visible in high layers casts 
doubts on their convective origin. Overall, we do not find indications of downflows that could be associated 
with overturning convection at our detection limit of 150 m s" 1 . Either no downflows exist, or we have been 
unable to observe them because they occur beneath r = 1 or the spatial resolution/height resolution of the 
measurements is still insufficient. 

Subject headings: convection - sunspots - Sun: photosphere - Sun: surface magnetism 
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1. INTRODUCTION 

The origin of the Evershed flow is not completely un- 



dersto o d dBorrerol I2009t iScharmeri I2009t iSchliche nmaier 
120091: iBellot Rubiol 120101: iNordlund & Scharmerl 1201 Oh . 
Models based on moving flux tubes explain the flow 
in terms of a pressure gradient that builds up along 
the tubes as they rise from the sunspot ma gnetopause 
(Schlichenmaie r et alj[l99 8: Schlic henmaieri 120021) . whereas 
siphon flow models invoke pressure gradients created by 
different field str engths at the footpoints of elevated mag- 
netic arches (e.g.jMever & Schmidi ll 9681: iDegenhardtll 19891 : 
iMontesinos & Thomasll 19971) . In both cases, mass conserva- 
tion is secured by downflows that occur in the mid and outer 
penumbra as the field lines return to the solar surface. 

Another possibility is that the Evershed flow is caused by 
convection in the presence of the sunspot magnetic field, 
which induces an aniso tropy in the radial di rection. The 
idea was put forward by Scha rmer et al.l d2008l) and seems to 
be supported by radiative magnetohydrodynam ic simulations 
dHeinemann et all 120071: iRempel et alj|2009allbl) . According 
to the simulations, hot weakly magnetized material ascends in 
the penumbra and becomes nearly horizontal after being de- 
flected outward by the inclined sunspot field. This results in a 
penumbral filament — an elongated overturning flow pattern 
with an upward component at the center (the Evershed flow) 
and lateral downflows at the edges. The process is similar to 
granular convection in the quiet Sun, except for the existence 
of a preferred horizontal direction. 

In the simulations one observes narrow (0"2-0"4) lanes of 
downflows on either side of the filaments. They appear where 



the Evershed flow returns back to the solar surface. The down- 
flows ensure mass conservation and reach velocities of up to 
1.5 km s" 1 at optica l depth r = 0.1 dHeinemann et al J 120071 : 
Remp el et alj |2009a). Interestingly, they have a small hor- 
izontal component toward the umbra, i.e., their direction is 
opposite to that of the Evershed flow. 

A clear detection of downward motions at the edges of 
penumbral filaments would support the existence of overturn- 
ing convection in sunspots. However, this is not an easy task 
because the Evershed flow is much stronger and may hide 
them in an efficient way, especially at low angular resolution. 
To minimize the problem, it is convenient to investigate the 
sunspot regions perpendicular to the line of symmetry, i.e., 
the line connecting the sunspot center with the disk center; 
there, the Doppler shifts induced by radial flows are zero, fa- 
cilitating the detection of vertical motions. For an unambigu- 
ous identification of downflows using Doppler measurements, 
the spot should be as close to the disk center as possible, so 
that redshifts can reliably be associated with downflows and 
blueshifts with up flows. 

Despite these problems, there have been some reports 
of small-scale convect i ve mo tions in sunspot penumbrae. 
Sanchez Almei da et al.l d2007l) found a positive correla- 
tion between upward velocities and brightness from high- 
resolution spectroscopic measurements taken at the Swedish 
lm Solar Telescope (SST). The observed correlation is remi- 
niscent of normal convection, suggesting that a similar mech- 
anism might be at work in the penumbra. Rim meld d2008l) 
studied penumbral flows near the disk center using the Uni- 
versal Birefringent Filter at the Dunn Solar Telescope. He de- 



2 



BELLOT RUBIO ET AL. 



tected upflows in the dark core and downflows on either side 
of a penumbral filament, as expected from overturning con- 
vection. However, other filaments in the same filtergrams did 
not show downflows. 

Using the Hinode satellite, Ichimoto et al. (2007) discov- 
ered an apparent twisting motion of brightness fluctuations in 
penumbral filaments located at ±90° from the symmetry line. 
The direction of the twist, as well as the associated Doppler 
signal (albeit at lower angular resolution), were consistent 
with transverse motions from the center of the filament to the 
edge facin g the observer. Si milar conclusions have been ob- 
tained by Bharti et al] (120101) from a larger sample of Hinode 
filtergrams. These authors measured twist velocities of more 
than 2.1 km s" 1 in the portion of the filaments closer to the 
umbra and about 1 km s" 1 at larger distances. IZakharov et al.l 
(2008) reported the detection of Doppler signals compatible 
with horizontal overturning motions in a filament observed 
almost perpendicularly to the line of symmetry, 40° away 
from the disk center. Their measurements were taken with the 
SOUP magnetograph at the SST The line-of-sight velocity of 
the horizontal motions was determined to be around 1 km s" 1 
from a Milne-Eddington inversion of the data. Unfortunately, 
the Fe I 630.25 nm line scans performed by Zakharov et al. 
required 123 s to be completed, meaning that any change in 
the flow field, solar scenery, or even seeing conditions dur- 
ing the scan could have affected the measurements. Velocities 
of 1 km s" 1 extending over 0"4 seem within easy reach and 
should have b een detected earlier. 

By contrast. [Franz & Sch lichenmaier (2009) did not find in- 
dications of overturning downflows in Doppler velocity maps 
computed from spectropolarimetric observations taken with 
Hinode at a resolu t ion of 0."3. Also the theoretical work of 
iBorrero & SolankH J2010) suggests that small-scale overturn- 
ing motions are not needed to explain the net circular polariza- 
tion observed in the penumbra, since the much stronger E ver- 
shed flow already accounts for it (see IBorrero et al. 2007). 

In view of these results, the existence of overturning con- 
vection cannot be considered as established. Thus, it is impor- 
tant to continue the search with measurements of the highest 
quality. Here we present an analysis of sunspot observations 
made with the Littrow spectrograph of the SST. The spot was 
located very close to the disk center, thus satisfying one of 
the requirements for good sensitivity to downward motions 
mentioned above. These observ ations reach a spatial reso lu- 
tion of about 0"2 and were used Bellot Rubio et al. (2005) to 
study the flow field of dark-cored penumbral filaments. To 
our knowledge, they are the highest spatial-resolution spec- 
troscopic measurements of sunspot penumbrae ever obtained. 
The resolution is sufficient to distinguish the dark cores of 
penumbral filaments, which are only barely detectable in the 
seeing-free observations of the Hinode spectropolarimeter. 



2. OBSERVATIONS AND DATA ANALYSIS 

NOAA Active Region 10 756 was observed with the 
Swedish 1 m Solar Telescope (Sch armer et al.ll2003b on 2005 
May 1. The best seeing conditions occurred at around 10:22 
UT, when the heliocentric angle of the spot was 5°. Figure Q] 
shows a full disk MDI continuum image taken at 10:17 UT, 
while Figure |2] displays a context G-band continuum filter- 
gram proc essed with the Multi-F rame Blind Deconvolution 
technique (Ivan Noort et al.l 120051) . The resolution of the G- 
band continuum image approaches 0"1. 




Figure 1. Full disk continuum image taken at 10:17 UT on 2005 May 1 by 
the Michelson Doppler Imager aboard SOHO. The spot closer to disk center 
is NOAA AR 10756, at coordinates S6.9° and W4.3 . Courtesy Stanford 
Lockheed Institute for Solar Physics. 

2.1. Spectroscopic observations 

We used the TRI-Port Polarimetric Echelle-Littrow (TRIP- 
PEL) spectrograph and the adaptive optics system of the SST 
to obtain line spectra at the highest angular resolution possi- 
ble. TRIPPELQ has a grating of 79 grooves mm a blaze an- 
gle of 63.43°, and a theoretical resolving power of 240000. 
Fe I 557.6 nm, Fe II 614.9 nm, and Fe I 709.0 nm were 
recorded simultaneously during the observations, but in this 
paper we consider only the 709.0 nm measurements because 
their short exposure times provided the best angular resolu- 
tion. 

Fe I 709.0 nm is a non-magnetic line suitable for Doppler 
shift measurements. Its na rrow shape makes it very se nsitive 
to line-of-sight velocities (Cabrer a Solana et al.l 120051) . We 
performed spectroscopic observations of the region around 
709 nm with a wavelength sampling of 1.046 pm, a pixel 
size of 0."04, and an exposure time of 200 ms. The spec- 
tra were recorded on a Kodak Megaplus 1.6i camera. The 
spectrograph slit had a width of 0"1 1 (25 /im) and covered a 
length of 37"8 (945 pixels). Simultaneously, slit-jaw images 
were taken through a 1 nm wide filter centered at 694 nm, also 
with exposure times of 200 ms. 

The main advantage of slit spectrographs over tunable fil- 
ters is that they record the full line profile at once, thus pre- 
serving spectral integrity. The disadvantage is that the solar 
image has to be stepped across the spectrograph slit to create 

1 See http://dubshen.astro.su.se/wiki/index.php/TRIPPEL_spectrograph 
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Figure 2. Left: G-band continuum image of AR 10756 obtained on 2005 May 1 at 10:22 UT through a 1.1 nm wide filter centered at 436.4 nm. The image has 
been reconstructed from five individual filtergrams of 20 ms exposure time. The arrow marks the direction to disk center (DC). The white line represents the 
spectrograph slit. It samples the umbra, the inner penumbra perpendicular to the symmetry line, and part of the center-side penumbra. Right: Slit-jaw image 
at 694 nm corresponding to the spectra discussed in this paper. The small dashes A and B indicate the portion of the slit shown in Figure[3] Two dark-cored 
penumbral filaments (PF 1 and PF 2) are clearly seen at the positions indicated by the arrows. 



2D maps. We performed 20"-wide scans of the spot, but only 
a few slit positions attained a spatial resolution better than 
0"25. For this reason we concentrate on the best slit position, 
not on the entire map. As can be seen in Figure 2, the slit 
crossed the center-side penumbra and the region perpendicu- 
lar to the line of symmetry near the umbra, providing a wide 
range of positions where overturning convection could occur. 

The spectral images have been corrected for dark current, 
flatfield, and optical distortions (smile and keystone). The lat- 
ter produce curvature of the spectral lines and a variation of 
the disp ersion along the sl i t, res pectively. More details can be 
found in Langange n et al.l (12007b . 

2.2. Velocity measurements 

We derive absolute velocities from the observed spectra as 
follows. First, we set up a relative wavelength scale using 
the dispersion obtained from the pixel separation between two 
lines with well-determined solar wavelengths in the corrected 
flatfield images. To calibrate this scale in absolute terms, it is 
necessary to fix the wavelength of one pixel in the spectrum. 
We use the telluric H2O line at 709.4050 nm for that purpose. 



The H2O line is present in the same spectral region and is not 
affected by motions in the solar atmosphere, so it provides an 
excellent reference. Its central wavelength has been obtained 
from th e Fourier Transform Sp ectrometer atlas of the quiet 
Sun by iBrault & Neckell (119871) with a precision better than 
50 m s -1 . For each height along the slit, we determine the 
position of the H2O line core by means of a parabolic fit. The 
derived values show pixel-to-pixel fluctuations of 110 m s -1 
due to the weakness of the line, but we reduce them down to 
40 m s" 1 applying a 12-pixel boxcar average along the slit. 
The final result is one position in each spectrum with a very 
precise wavelength (that of the H2O line), i.e., an absolute 
wavelength scale. 

The spectrally resolved profiles of Fe 1 709.0 nm are used to 
derive Doppler shifts by means of line bisectors. We compute 
bisectors for intensity levels from 0% to 88%, where 0% rep- 
resents the line core and 100% the continuurrQ, through linear 
interpolation of the original wavelength samples. The bisec- 

2 We do not consider intensity levels higher than 88% to avoid the CN 
blend present in the very far red wing of the Fe I line (at 709.069 nm). 
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Figure 3. Absolute line-of-sight velocities derived from the line bisectors at different intensity levels: 0-8% (thick black line), 50-58% (thin black line), and 
80-88% (red line). Positive velocities indicate redshifts. Since the spot was very near the disk center, they most likely correspond to downflows. The horizontal 
solid and dashed lines mark velocities of and 150 m s , respectively. The dotted blue line displays the continuum intensity along the slit, in a relative scale. 
The vertical lines labeled PF1 and PF2 in the upper abcissa mark the position of the two dark-cored penumbral filaments crossed by the slit (see Figure|5J. 



tor positions are converted into Doppler shifts by subtracting 
the laboratory wavelength of the Fe I line (709.03835 nm ac- 
cording to Nave et al. 1994), and then transformed into line- 
of-sight velocities. 

Different intensity levels sample different layers of the at- 
mosphere. Generally speaking, higher bisector intensities cor- 
respond to deeper layers . In the penumbral atmosphere of 
iBellot Rubio et al. (2006), the far line wing (as represented 
by intensities between 80% and 88%) is formed in the first 
100 km above optical depth unity, i.e., the bisector Doppler 
shifts measured in this range reflect flow velocities from the 
very deep photosphere. 

The bisector velocities still need to be corrected for gravi- 
tational redshift (636 m s" 1 ) and relative motions between the 
Sun and the observer . The former, computed as explained by 
iMartmez Pillet et alJ (119971) . amount to +316 m s" 1 , of which 
-261 m s" 1 correspond to the Earth's rotation, +436 m s" 1 to 
the Earth's orbital motion, and +141 m s" 1 to the solar rotation 
at the position of the spot (negative velocities are blueshifts). 

Our absolute velocity calibration does not rely on the gran- 
ulation or the umbra, and therefore it is not affected by uncer- 
tainties in the convective blueshift of the line or the presence 
of umbral flows. The main source of systematic errors comes 
from the laboratory wavelength of Fe I 709.0 nm , which has 
an uncertainty of about 100 m s" 1 dNave et alJl994h . Another 
source of systematic error is the central wavelength of the ref- 
erence H2O line (50 m s -1 ). These errors add linearly, so the 
maximum systematic error of our calibration is 150 m s" 1 . 
The effect of a systematic error is to shift all the velocities 
up or down as a whole. In contrast, the pixel-to-pixel "noise" 
of the velocity curves is due to random errors. We estimate 
the random error to be on the order of 1 10 m s" 1 . This error 
comes from uncertainties in the position of the reference H2O 
line (40 m s" 1 ) and uncertainties in the determination of the 
bisector position (about 100 m s" 1 ), added quadratically. All 
in all, our velocity measurements should be accurate to within 
±110ms _1 , with systematic errors below 150 ms" 1 . 

3. RESULTS 

Figure [3] show the Doppler velocities observed along the 
slit, excluding the dark umbra where it is not possible to com- 



pute reliable positions for the weak H2O line. The outer part 
of the center-side penumbra is to the right. Displayed are 
line-core and line-wing velocities corresponding to the bisec- 
tor shifts averaged between the 0% and 8% intensity levels 
(thick black line), between 50% and 58% (thin black line), 
and between 80% and 88% (red line). As mentioned before, 
higher intensity levels progressively sample deeper layers of 
the photosphere. For reference, the horizontal lines indicate 
velocities of and 150 m s" 1 (positive values represent red- 
shifts). The continuum intensity is shown in the top part of 
the figure to help identify bright and dark structures (dotted 
blue line). 

The first thing to note from Figure [3] is the tendency of all 
the velocities to increase toward more blueshifted values as 
the outer penumbra is approached, i.e., from left to right. This 
is due to two reasons: the strong enhancement of the magni- 
tude of the Evershed flow with radial distance from the spot 
center (which compensates the increasing inclination of the 
flow), and the more favorable projection of the flow velocity 
to the line of sight toward the upper end of the spectrograph 
slit. Indeed, the angle between the penumbral filaments and 
the line of sight decreases from about 90° near the umbra to 
about 40° at the top of the slit; since the Evershed flow occurs 
along the filaments, the projection leads to stronger Doppler 
velocities. At an heliocentric angle of only 5°, however, the 
first effect is dominant. 

We also note that the bisector velocities corresponding to 
intensity levels of 0%-8% and 50%-58% are very similar. 
This means that the lower half of the line show essentially 
vertical bisectors (i.e., no line asymmetries). By contrast, the 
velocities are substantially larger and to the blue in the far 
wings. Because of the larger velocities at high intensity lev- 
els, the bisectors are tilted to the blue near the continuum. Fig- 
ure |4j a) displays a typical example: the strong tilt is caused 
by a "satellite" that starts to be visible in the far blue wing of 
the line. Also depicted in the figure are bisectors showing a 
more gradual (but strong) shift to the blue near the continuum, 
as well as cases of vertical bisectors with zero velocities and 
bisectors showing a tilt to the red. The latter, however, are not 
common in this part of the penumbra. 

Deep in the atmosphere, as sampled by the bisectors near 




Figure 4. Examples of Fe I 709.0 nm line profiles (top) and bisectors (bottom) at different positions along the slit, (a) Vertical bisectors with a strong tilt to the 
blue near the continuum, (b) Strongly tilted bisectors showing larger blueshifts toward the continuum, (c) Vertical bisectors with zero velocities, (d) Bisectors 
tilted to the red toward the continuum. These examples correspond to pixel positions 338, 54, 209, and 552 in Figure 3, respectively. 



the continuum, we always find blueshifts or zero velocities 
along the slit. The maximum blueshifts reach 1 .2 km s" 1 and 
tend to be cospatial with bright penumbral filaments, although 
we sometimes see a small displacement of the maximum ve- 
locity away from the observer, relative to the peak brightness 
position. In the dark regions outside of the filaments we also 
detect blueshifts for the most part. However, the velocities 
are strongly reduced. Sometimes they even drop to zero (e.g., 
Figure|4j;). Although these instances are not common, this is 
the first time that zero velocities are observed in dark areas of 
the center-side penumbra. Earlier measu rements at lower spa- 
tial resolution did not reveal them ("e.g.. iHirzberger & Kneerl 
l200lt iRouppe van der Voortl l2002; Bell ot Rubio etail l2006). 
The blueshifts depicted in Figure[3]may represent radial flows 
along inclined flux tubes. In that case, one would expect 
very small or no plasma motions in between the flow chan- 
nels. However, this is not observed. The relative absence 
of zero velocities in our data may indicate that the tubes are 
still partly unresolved, or that there is some amount of stray 
light contamination. Another possibility is that the volume 
outside the flow channels is not c ompletely at rest, as sug- 
gested by Stokes i nversions (e.g., iBellot Rubio etail 12001 
iBorrero et alj 120051) and net circular po larization measure- 
ments with Hinode dlchimoto et alj|2008l) . 

The line-core shifts (thick black curve in Figure[3]l are much 
smaller than their line-wing counterparts. They do not exceed 
200 m s" 1 except in the middle center-side penumbra. Most of 
the positions along the slit show blueshifts. Redshifts do also 
occur, but with velocities of less than 100-150 m s" 1 . These 
values are smaller than the systematic errors of our absolute 
velocity calibration. Apparently, the redshifts do not bear any 
relationship to the penumbral filaments: sometimes they oc- 
cur on one side of the filament, sometimes on the other. The 
only exception is perhaps the structure located at pixel po- 
sition 20, which shows relatively constant redshifts of about 
100 m s" 1 . This filament makes an angle of about 90° to the 
symmetry line. 



It is also important to note that the dark-cored penumbral 
filaments crossed by the slit (PF 1 and PF 2 in Figure [3} do 
not exhibit blueshifts in the dark lanes and redshifts in the two 
lateral brightenings, contrary to what would be expected from 
a resolved overturning flow structure. 

4. DISCUSSION 

The observations analyzed in this paper are very favor- 
able for the detection of downflows in the penumbra because 
they correspond to a sunspot located only 5° away from disk 
center. The contribution of horizontal flows to the observed 
Doppler signals is small, and redshifts/blueshifts can safely be 
associated with vertical downflows/upflows. Moreover, our 
observations provide full line profiles at the highest angular 
resolution achieved in this kind of measurements (better than 
0."25), including a telluric line that we have employed to cal- 
ibrate the velocity scale in absolute terms. 

Despite the excellent quality of the data set, we do not de- 
tect downflows that could be associated with overturning con- 
vection in deep layers: near the continuum the measured ve- 
locities are always to the blue. Higher in the atmosphere, as 
sampled by the line core, we observe some regions that could 
harbor downflows, but with small velocities not exceeding 
150 m s" 1 . The fact that they are visible only in high lay- 
ers suggests that they are not related to convective processes 
in the deep photosphere. They could be produced by down- 
flows associated with the inverse Evershed flow, but also by 
penumbral oscillations or waves. In the absence of time series 
we cannot decide between the different scenarios. 

Overall, our data suggest that downflows due to overturning 
convection are not larger than 150 m s" 1 in the photospheric 
layers accessible to the observations, while the simulations 
predict values up to 1.5 km s" 1 . The lack of a clear detec- 
tion of overturning downflows may simply indicate that they 
do not exist. However, it is also possible that they have gone 
undetected if (a) they occur beneath r = 1 ; (b) they have veloc- 
ities smaller than 150 m s" 1 ; (c) they are concentrated in very 
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thin sheets not resolved by our observations; or (d) they fill 
only a small fraction of the formation region of the 709.0 nm 
line. 

To examine possibility (c) we need spectroscopy at a reso- 
lution of 0."1 or better. Even with infinite spatial resolution, 
(d) may hamper the detection of the downfiows if their vertical 
extent is much smaller than the width of the contribution func- 
tions of typical photospheric lines. In fact, the simulations 
show that the downflow lanes surrounding the penumbral fil- 
aments are not completely perpendicular to the solar surface, 
so that at some points they have a thickness o f only 50-100 
km in the vertical direction (e.g., Figure 9 of Remp el et alJ 
2009a). Such small structures may not be able to leave clear 
signatures in the emergent intensity profiles. Line synthesis 
calculations based on the simulation results are needed to clar- 
ify this aspect. If Fe I 709.0 nm does not provide sufficient 
height resolution, then the search for overturning downfiows 
should be continued using lines with narrower contribution 
functions, like C I 538.03 nm. Another requirement to make 
progress is to increase the realism of the simulations. Current 
simulations reproduce many of the observed characteristics of 
the penumbra, but they predict downfiows that we do not de- 
tect in spite of our excellent spatial, temporal, and spectral 
resolution. Thus, it is important to confirm that their results 
are not affected by too large values of the viscosity and mag- 
netic diffusivity, by inadequate boundary conditions, or by too 
shallow computational boxes. Hopefully, the combination of 
these observational and numerical efforts will lead to a better 
understanding of the penumbra and the Evershed flow. 
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